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Abstract 

We study the CMB constraint on non-Gaussianity in CDM isocurvature pertur- 
bations. Non-Gaussian isocurvature perturbations can be produced in various mod- 
els at the very early stage of the Universe. Since the isocurvature perturbations little 
affect the structure formation at late times, CMB is the best probe of isocurvature 
non-Gaussianity at least in the near future. In this paper, we focus on uncorrelated 
isocurvature perturbations and constrain their non-Gaussianity. For this purpose, 
we employ several state-of-art techniques for the analysis of CMB data and simula- 
tion. We use the WMAP 7 year data of temperature anisotropy. When the adiabatic 
perturbations are assumed to be Gaussian, we obtained a constraint on the isocur- 
vature non-Gaussianity a 2 f^°^ = 40 ±66 for the scale invariant isocurvature power 
spectrum, where a is the ratio of the power spectrum of isocurvature perturbations 
to that of the adiabatic ones. When we assume that the adiabatic perturbations can 
also be non-Gaussian, we obtain /nl = 38 ± 24 and a 2 f^°^ = —8 ± 72. We also 
discuss implications our results on the axion CDM isocurvature model. 



1 Introduction 



The adiabaticity, or isocurvature mode of primordial density fluctuations is one of impor- 
tant probes of cosmology in various respects. Although current cosmological observations 
suggest that primordial density fluctuations are almost adiabatic and they give severe 
constraints on the size of isocurvature fluctuations, some fraction of their contribution 
is still allowed [1]. Residual isocurvature fluctuations can be generated when there exist 
multiple components with different origins, which are associated with dark matter, baryon 
and neutrino [2] #1 . Such examples include axion [5, 6, 7, 8, 9, 10, 11, 12] and Affleck- 
Dine baryogenesis [13, 14, 15], where cold dark matter (CDM) and baryon isocurvature 
modes can be respectively generated. These isocurvature modes are basically uncorre- 
lated with adiabatic ones, however, when one considers a scenario where a light scalar 
field other than the inflaton is responsible for (adiabatic) density fluctuation such as the 
curvaton model [16, 17, 18], isocurvature perturbations can be correlated with the adia- 
batic ones and be easily generated, depending on how and when CDM and baryon are 
created [19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31]. In addition to above men- 
tioned examples, a variety of models with isocurvature fluctuations has been extensively 
studied in the literature, hence the information on isocurvature fluctuations would give a 
lot of insight on various aspects of cosmology, particularly on models of dark matter and 
baryogenesis as well as those of the early Universe. 

Although the theoretical works on non-Gaussianity from isocurvature fluctuations have 
been relatively well studied [32, 33, 34, 35, 36, 37], observational constraints on them have 
not been explored much* 2 . In particular, we cannot find any work investigating this issue 
by using the actual data except Ref. [48], where the constraint has been studied with 
Minkowski functionals using the WMAP 3 data, although there are a few papers which 
elaborate its future CMB constraints [49, 50, 51]. In the light that we now have precise 
cosmological data to explore non-Gaussianity as seen from the counterpart for adiabatic 
ones, non-Gaussianity in isocurvature fluctuations should be pursued more. 

In this paper, we investigate the optimal constraint on non-Gaussianity in CDM 
(baryon) isocurvature fluctuations from CMB bispectrum estimator. Although, as men- 
tioned above, isocurvature fluctuations can be correlated with adiabatic ones in some cases 
and there are also other kinds of modes such as neutrino one, we in this paper present the 
methodology of our analysis and concentrate to report the constraint on the CDM (baryon) 
uncorrelated type. Constraints on other types such as correlated ones and neutrino modes 
will be reported in a forthcoming paper [52]. 

The organization of this paper is as follows. In the next section, we give the for- 
malism to discuss non-Gaussianity in models with isocurvature fluctuations and also set 
our notation. In Section 3, we describe our analysis method to obtain the constraint on 
non-Gaussianity from isocurvature fluctuations. Then in Section 4, we present our results 

* 1 Isocurvature fluctuations in dark radiation have also been studied in Refs. [3, 4]. 
# 2 At linear order, there are a lot of works which study the observational constraints [38, 39, 40, 41, 42, 
43, 44, 45, 46, 47]. 
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on the constraint. As an application of our results, we consider constraints on the axion 
isocurvature model in Section 5. The final section is devoted to the conclusion of this 
paper. 

2 Model of non-Gaussian perturbations and CMB 
signature 

We consider primordial curvature and CDM isocurvature perturbations, respectively de- 
noted as ( and S, in the following form: 

m = cg(x) + jjwcco?) 2 - (Cg(^) 2 », (i) 

S(x) = S G (x) + fSl°\S G (x) 2 - (S G (x) 2 )), (2) 

where ( G and S G are Gaussian parts of the primordial curvature and isocurvature per- 
turbations, respectively, /nl and f^i°^ are the non-linearity parameters of the curvature 
and isocurvature perturbations, respectively. In the following, we denote these primordial 
perturbations with X A (x). Then Eqs. (1) and (2) can be recast into 

X A (x) = X A (x) + fUXoW 2 - (X£(x) 2 }), (3) 

where X A is the Gaussian part of X A . Here we defined a non-linearity parameter f$ L , 
which is related to the adiabatic and isocurvature ones via f^ L = |/nl and f^ L = f^i°^- 
In this paper, we consider uncorrelated isocurvature perturbations, so that (( g Sg) — 0. 
Thus only correlation functions which contain either ( or S are non-zero. In terms of 
the Fourier components of X A , the bispectrum from either the primordial curvature or 
isocurvature perturbations is 

(X^X^X^h)} = 2f A L (P X A(h)P xi (k 2 ) + (2 perms)), (4) 

where P X A{k) is the power spectrum of the Gaussian perturbations X A . Here we neglected 
loop contributions although we take those into account in Sec. 5. 

Neglecting the secondary non-Gaussianities arising from the second or higher order 
cosmological perturbation theory, the harmonic coefficients of the CMB temperature 
anisotropy from primordial perturbations X A can be given as 

j^gf(k)X A (k)Y lm (k), (5) 

where gf(k) is the temperature transfer function for the primordial perturbations X A . 
The total CMB anisotropy is the sum of those from the curvature and isocurvature per- 
turbations, i.e. ai m = a\ m + af m . 
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Since ( and S are uncorrelated we need to consider the polyspectra which contain only 
either ( or S. Then the total polyspectra are the sum of those from each perturbation. 
Angular power spectrum Cf, which is defined by (atm a Vrn*) = Q^H'^mm' can be given as 

Cf = lj dkk 2 gf(k) 2 P x ,(k). (6) 

The reduced bispectrum &; 1 i 2 / 3 is defined by 

(7) 

where ^* m3 is the Gaunt integral, which can be written in terms of the Wigner-3j 
symbol as 



h i 2 i 3 _ ,/ (2/ 1 + l)(2/ 2 + l) (2/ 3 + l) / h h k \ (h h h 



G 

y mi m 2 m 3 y 4?r m2 ms j ^ooo 

Given a bispectrum of primordial perturbations in Eq. (4), reduced CMB bispectrum from 
each of the primordial perturbations can be given as 

hf ll2l3 = 2f4 L J drr 2 « (r)^(r)^(r) + (2 perms)) , (9) 

where af (r) and /3f(r) are 

= lj dkk 2 gf(k) 3l (kr) } (10) 
tf(r) = IJ dkk 2 gf(k) 3l (kr)P X A(k). (11) 
For later convenience, we introduce a normalized bispectrum 

The amplitude of isocurvature power spectrum is constrained from CMB angular power 
spectrum. In Ref. [1], it is shown that the WMAP observation of CMB gives an upper 
bound on the fraction of isocurvature power spectrum in the total one 

a = ^<0.15 (13) 

at 95 % confidence level. 

Note that the bispectrum is proportional to a 2 f^°\ This can be seen from Eqs. (9) 
and (11). Hence we report our constraint on non-Gaussianity in isocurvature fluctuations 
for the combination of a 2 f^°\ 
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3 Analysis method 



Here we describe our analysis method to give the constraint on the non-linearity parameters 
from CMB data. Our method basically follows the one in Ref. [54]. 



3.1 Estimator of non-linearity parameters 

In the limit of small non-linearity parameter f£ L , the effect of the deviation from Gaussian- 
ity manifests in the CMB bispectrum, so that the optimal estimator can be constructed 
from the three point function of CMB anisotropy [55]. We adopt the following cubic 
estimator of non-linearity parameters 



/NL / \ ( ^prim )fL=l ( ^prim ) /| L = 




(14) 



where the square bracket with subscript f£ L = 1 indicates an ensemble average over 
simulations with non-zero non-linearity parameter shown in the subscript with the other 
non-linearity parameter being fixed to zero. S^ rira is the cubic statistics [55, 56, 57, 58], 
which is given by 

)o] , (15) 

{lm} 

where the square bracket with subscript indicates an ensemble average over Gaussian 
simulations. a\ m is a harmonic coefficient obtained from observed (or simulated) data 
maps with suitable weighting, which will be discussed in Section 3.3. Eq. (14) can be 

schematically represented as f$ L = y^(5'^. im )~i =1 Sp Tim . 

B NL 

In our analysis, to estimate {ai m aii m >) we accumulated 500 Monte Carlo (MC) samples. 
This number is sufficient, since the uncertainty in estimation of due to this limitation 

of MC samples is much smaller than the Gaussian error yj (Sp TiTa Sp Vim ) . 



Assuming the Gaussianity for X A , we estimate the covariance of the estimator ff- 



NL; 



(/nl/nl)o — (Sprim) ja' =] (^rim'^prim)o (^prim) fB' =] 

(16) 

A'B' NL ' NL 

(^prim) fB 



Here we used the relation (Sp rim ) jb =1 = (Sp rim Sp Tim )o- See Appendix A for the derivation. 

The estimator in Eq. (14) can be regarded as the generalization of the fast estimator in 
Ref. [55] to the case where the initial perturbations are mixture of non-Gaussian adiabatic 
and uncorrelated isocurvature perturbations. When we assume that either of adiabatic 
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or uncorrelated isocurvature perturbations are non-Gaussian and the other is Gaussian, 
Eq. (14) can be reduced to 

/nL = ^prim/ (^prim) /n l =1> K^J 

where the subscripts A indicates the perturbations which are assumed to be non-Gaussian. 
In this case, the variance of f$ L is given by 1/ (S^ Tim ) 1. 

We note that there is a difficulty in computing the normalization factor (S^. im ) fB =1 . If 
we naively evaluate Eq. (15) with small but non-zero f£ L and average them, a substantial 
number of simulations are required due to a large Gaussian fluctuation. Instead, we divide 
di m into its Gaussian and non-Gaussian parts and evaluate «S^. im without terms which are 
to vanish by averaging. The details are presented in Appendix B. This treatment also 
removes the need for setting a non-zero fiducial a. 



3.2 Non- Gaussian CMB simulation 

In order to determine the normalization (S , ^ rim ) f s _i, we need to simulate non-Gaussian 

\ pi mi/ Jp^L — 

CMB maps. With the Fourier transformation, Eq. (5) can be rewritten as 

<*L = J drr 2 a?(r) J dfY lm (r)X A (r). (18) 

In the case of local-type non-Gaussianity, Eq. (18) allows us to simulate the non-Gaussian 
CMB maps exactly [59, 60]. We adopt the fast method developed in Ref. [61], which 
enables simulations of the local-type non-Gaussianity with sufficient speed. The key tech- 
nique in the method of Ref. [61] is that the line of sight integral in Eq. (18) is evaluated 
by the Gaussian quadrature with optimized nodes and weights. In our simulation, we 
optimized the nodes so that the mean square of the error in simulated maps a\ m at each 
multipole (Im) should be less than 0.01. For Z max = 1024, we found that this level of accu- 
racy requires 42 and 15 nodes for curvature and isocurvature perturbations, respectively. 



3.3 Optimally-weighted CMB maps 

In estimating f^ L , we need to suitably weight observed (and simulated) maps ai m to obtain 
di m , m order to take into account the sensitivity and resolution of the survey. As discussed 
in Ref. [54] , the optimal weight is the inverse of the variance. This can be schematically 
represented as di m = [(7 -1 a]i m , where C = C$ + Cjq is the covariance matrix, with those 
of signal and noise being denoted as Cs and CV, respectively. While Cs is diagonal in the 
harmonic space, Cn is in general complicated for real observations. In particular, WMAP 
and many CMB surveys have multiple channels, so that we need to take into account 
different beam functions, inhomogeneous noises, and a sky cut. In such a case, d\ m can be 
given in an implicit form as 

(Cs 1 + C N l )C s d = C N \ (19) 
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where = b^Cff and C Ar 1 a = cffi a^, with the subscript z indicating a 

i i 

channel of the survey. feW here is the beam function. 

Since the matrix of (C^ 1 + C^ 1 ) is dense in both the harmonic and pixel spaces, 
direct implementation of inversion (Cg 1 + C^ 1 ) -1 is substantially prohibited. Instead, 
Eq. (19) can be solved by the conjugate gradient (CG) method with preconditioning 
[62]. How fast a CG method converges crucially depends on a choice of a pre-conditioner. 
We adopt a fast method developed in Ref. [63], which makes use of the multi-grid pre- 
conditioning. This method also allows to marginalize over the amplitude of components 
which are proportional to template maps {r}. To do this, C^ 1 is replaced with C^ 1 — 

C] v 1 T a (T a C N 1 Ti ) )~ 1 C N 1 Ti J , where the subscripts a and b indicate template components. 

ab 

For the details of the method, we refer to Refs. [63, 54]. 

4 Result 

In our analysis, we assume a flat power-law ACDM model as a fiducial model and adopt 
the mean values for cosmological parameters from the WMAP 7 year data alone [1]. We 
also assume that the power spectrum of isocurvature perturbations is in a power-law 
form Ps G (k) oc k n ' ls °~ i 1 and with regard to the fiducial value of the spectral index n; so we 
consider two different C3jSGS^ ^'iso 

= 0.963(= n a di) and n iso = 1. The transfer function of 
CMB is computed using the CAMB code [64]. We combine the foreground-cleaned maps 
of V and W bands of the WMAP 7 year data [65, 66] #3 with a resolution A^ side = 512 
in the HEALPix pixelization scheme [67] * 4 . We adopt the KQ75y7 mask [66], which 
cuts 28.4 % of the sky. We also set the maximum multipole Z max to 1024 in our analysis. 
We marginalize the amplitudes of the monopole I = and dipoles I = 1 as default. We 
also optionally marginalize the amplitudes of the Galactic foreground components at large 
angular scales using the synchrotron, free-free and dust emission templates from Ref. [66]. 
Our results are summarized in Table 1. 

To check our analysis method, we have evaluated a constraint on the non-linearity 
parameter for adiabatic perturbations assuming isocurvature perturbations being absent, 
and compare it with those in a previous study. In this case, we obtain /nl = 31 ± 21 at 
la level without template marginalization of the Galactic foregrounds. Without template 
marginalization, our central value by about 0.5<r deviates from that of Ref. [1], which gives 
/nl — 42 ±21 from the foreground cleaned V+W band data with a resolution N side = 1024. 
Since there are substantial differences in our analysis and that of Ref. [1], including 
different fiducial cosmological models and resolutions of the maps used, we believe that 
this level of difference is acceptable and our result is consistent with the previous study. 
With template marginalization, we obtain /nl = 32 ±21, which is exactly the same as the 

# 3 http://lambda.gsfc. nasa.gov 
# 4 http: / /healpix.jpl. nasa.gov 
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Setups 


t 

/NL 


n 2 f (ISO) 


n iso = 0.963 w/o template marginalization 


31 ± 21 


5 ± 63 




(36 ± 23) 


(-39 ±69) 


w/ template marginalization 


32 ±21 


40 ±63 




(32 ± 23) 


(0 ± 70) 


?T,i SO = 1 w/o template marginalization 


31 ±21 


19 ±62 




(34 ± 23) 


(-22 ± 70) 


w/ template marginalization 


32 ±21 


40 ±66 




(38 ± 24) 


(-8 ± 72) 



Table 1: Constraints on / NL and « 2 / NL for different setups. We adopted four different 
setups regarding the fiducial value of rii SO and template marginalization of the Galactic fore- 
grounds. Constraints without parenthesis are estimated by fixing the other non-linearity 
parameter to zero. On the other hand, ones with parenthesis are estimated by marginal- 
izing over the other non-linearity parameter. 



one given by the WMAP group [1]. We found that template marginalization little affects 
constraints on /nl- 

Now we present constraints on non-Gaussianity in uncorrelated isocurvature pertur- 
bations. We first assume that the adiabatic perturbations are Gaussian and fix /nl to be 
zero. For the cases of n iso = 0.963 and n iso = 1, we respectively obtain a 2 f^°^ = 5 ± 63 
and 19 ±62 at la level without template marginalization. With template marginalization, 
these change to 40 ± 63 and 40 ± 66. We found that the constraints on a 2 f^l°^ are not 
strongly affected by the fiducial value of rii S0 . On the other hand, the constraints are more 
or less dependent on the treatment of the Galactic foregrounds. The central values can 
change by 0.5a while the error remains almost unchanged. However, this shows that the 
effects of residual foregrounds are not severe. 

When we assume that both adiabatic and isocurvature perturbations can be non- 
Gaussian, we give a joint constraint on (/nl,Q! 2 /nl )• F° r the cases of n iso = 0.963 and 

n iso = 1, we respectively obtain (/ NL , a 2 / N L 0) ) = ( 36 ± 23 > ~ 39 ± 69 ) and ( 34 ± 23 > ~ 22 ± 
70) without template marginalization. With template marginalization, these changes to 

(/nl, a 2 /NL 0) ) = ( 32 ± 23, ± 70) and (38 ± 24, -8 ± 72). The error for each non-linearity 
parameter here is estimated by marginalizing over the other non-linearity parameter. In 
Fig. 1, we show 2D constraints in the /nl-« 2 /nl°^ plane for the cases with template 
marginalization. Due to the correlation of a 2 f^°^ with /nl, a simultaneous fit for both 
of these variables changes the central value of a 2 f^°\ which are not so constrained as 
/nl- We found that the dependence on n iso is weak and the central values of ck 2 /nl can 
change by 0.5a by the treatment of the Galactic foregrounds. 

Having all these results, we conclude that CMB data is consistent with Gaussianity 
at 2a level, even if the uncorrelated CDM isocurvature perturbations are included in 
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Figure 1: Joint constraints on /nl and a 2 f^°^ at 1 and 2a levels with template marginal- 
ization. Solid red and dashed green lines respectively show the cases for n iso = 0.963 and 
ni SO = 1, with the central values being indicated by the red plus and green cross. Blue star 
corresponds to the case where the perturbations are purely Gaussian. 

consideration. 

Let us compare our results with other studies. While our method is optimal based on 
the bispectrum, constraints on the same type of non-Gaussianity is studied in Ref. [48] 
based on the Minkowski functionals, which gives a 2 f^°^ = — 15 ± 60, when the adiabatic 
perturbations are assumed to be Gaussian and n iso = 1. We found our constraint is 
consistent with the previous study. On the other hand, we cannot find any improvement 
in the constraint, although our method is based on the optimal estimator and should 
be better than suboptimal ones. This suggests that the Minkowski functional method 
is almost optimal for uncorrelated CDM isocurvature perturbations. This can also be 
expected from the Fisher matrix analysis in Ref. [49], which gives a Cramer- Rao bound 
A(a 2 /N L )=60. 

As a joint constraint on / NL and f^i° > our results are the first one obtained from 
observed data. On the other hand, the same constraint is forecasted using the Fisher 
matrix analysis in Ref. [49]. We note that la errors of our results are consistent with the 
forecast. 
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5 Application to the axion model 



In this section, we apply the constraints on the isocurvature non-Gaussianity to the axion 
isocurvature model. First, we shortly describe how non-Gaussian isocurvature perturba- 
tions in axion CDM arise in the inflationary Universe based on the 5A^-formalism [68, 69]. 
The baseline of our derivation is the same as in Ref. [32] (See also Ref. [48]). Throughout 
this section, we denote the energy density of a component % on the uniform density hyper- 
surface of the total matter with Pi(x). We adopt the non- linear definition of isocurvature 
perturbations in a component i, 

Si(x) = 3(6-0(2), (20) 

where Q and ( are respectively the curvature perturbations on the uniform density hyper- 
surfaces of the component i and the total matter. According to the 5iV-formalism, p(x) 
should be given by p%{x) = p^ 1 ^ 1 '^ 10 ', where pi is the mean of Pi(x) and Wi is the 
equation of state of the component i. 

The axion is a pseudo Nambu-Goldstone boson of the Peccei-Quinn (PQ) U(l) sym- 
metry, which solves the strong CP problem in quantum chromodynamics (QCD). If the 
PQ symmetry is broken during inflation, the axion has a classical field value a, = F a 9 
and a vacuum fluctuation Sa, where F a is the axion decay constant and 9 G [— tt, 7r] is the 
initial misalignment angle during inflation. At high temperature, the axion is massless. 
As temperature T decreases the QCD phase transition takes place. At this moment, the 
axion becomes massive and starts the coherent oscillation around the true vacuum. This 
oscillation of the axion contributes to the energy density of CDM. We assume that CDM 
in the Universe is a mixture of the axion and other CDM components which are adiabatic. 

The coherent oscillation of the axion synchronously starts on the uniform energy density 
hyper-surface of the total matter at around m axion (T) ~ 3H(T) # 5 , where m axion and H 
is the mass of axion and Hubble expansion rate, respectively. The energy density of 
the coherent oscillation is proportional to the square of its initial amplitude, p ax i on (x) cx 
(aj + 5a(x)) 2 , which leads 



e s„u»(2) = 1 + 2 a ^ Q ( f ) + ^a 2 (x) - (5a 2 ) ^ 
a 2 a 2 



where 

al = a 2 + (6a 2 ). (22) 

Since other CDM components are assumed to be adiabatic, their energy density is uniform 
on the uniform density hyper-surface of the total matter. Therefore Pcdm(2) should be 
given by 

Pcbm(x) = Pcdm [(1 - r) + re 5axion ] , (23) 



# 5 Here we assume that there are no isocurvature perturbations in neutrinos or, if any, extra radiations 
[3] . Otherwise, the coherent oscillation does not start synchronously on the uniform density hyper-surfacc 
and as a consequence, additional isocurvature perturbations in the axion can be induced. 
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where r = p ax ion/pcDM is the energy fraction of the axion in CDM and (1 — r) is that of 
other CDM components. Thus we finally obtain the isocurvature perturbations in CDM, 



CDM 



In 



PCDM 



In [l — r + re 5axion ] 



2r^- + 



4r a? 



2r- 



a,<5a 



2r- 



QiSa 



+ ■■■ (24) 



The last equality is approximately valid for r <C 1. In the following, we keep the terms up 
to the second order. 

Adopting the value of CDM abundance f2cDM^ 2 = 0.11 estimated by current observa- 
tions [1], r is given by [70, 32] 



r = U 



10 12 GeV 



-0.82 



10 12 GeV 



9 2 + 



#inf/27T 

10 12 GeV 



(25) 



where the first and second terms in the square bracket correspond to contributions from the 
classical field eij and the fluctuation 5a, respectively. The fluctuation of axion 5a is almost 
scale-invariant (See also Ref. [12]) and its root mean square is given by a/ (5a 2 ) = H- m f/2ir, 
where H m f is the Hubble scale during inflation. Then power spectrum and bispectrum of 
CDM isocurvature perturbations are given by* 6 



(SCDM^OSCDM^)) = 4r 



(H iBi /2n) 



(FJ) 2 + (iW2vr)2 k\ 



2tt 2 

-(2n)W 3 \k 1 + k 2 ), (26) 



(5 , cDM(^i)5'cDM(/i;2)5 , cDM(/i;3)) = 8r 3 

'(2tt 



X 



(F a #) 2 + (# inf /27r) 2 J 

2\2 



(27) 



k 3 Jc 3 



+ (2 perms) 



(27if5 i3 \k 1 + k 2 + k 3 



Having all these results, constraints should be 



U 
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47 



10 12 GeV 

1.64 



-0.82 



10 12 GeV 



9 2 + 



10 12 GeV 



10 12 GeV 



-2.46 



10 12 GeV 



10 12 GeV 



e 2 + 



gjnf/27T_ Y 

10 12 GeV J 



2 2 / H mi /2TT X " 

' 1 10 12 GeV 



gmf/27T_ Y 

10 12 GeV J 

H inf /2n 
10 12 GeV 



H ini /27i 



< 1, 



< 0.15A S , (28) 



< UOAi 



^ 6 Although, in the analysis in the previous sections, we omitted the loop contributions, we include 
them here. We also assume that logarithmic factors appear in the expression are considered to be 0(1) 
and we here set them to be unity. 
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where A s = 2.43 x 1CT 9 [1] is the amplitude of the adiabatic perturbations. From top to 
bottom, these three equations correspond to the constraints on the abundance r < 1, the 
isocurvature power spectrum a < 0.15 (95 % CL) [1], and non-Gaussianity \oc 2 f^i°\ < 140 
(2cr level) from our results. We should, however, note that since loop contributions are 
omitted in deriving the constraints in the previous sections, thus our results cannot in a 
rigorous manner be applied to cases for F a 9 < if in f/27r, where loop contributions dominate 
polyspectra. Therefore, in the parameter regions with F a 9 < H in f/2ir, our constraints 
should be deemed as rough estimate. 

In Fig. 2, we plotted the constraints of Eq. (28) in the H in f-9 plane with several fixed 
values of F a . Constraints have critical points at F a 9 = H in f/2ii, where contributions of the 
classical field value and the fluctuation are comparable. We found that the constraint on 
non-Gaussianity from our results gives an upper bound on H m { comparable to that from 
one on the isocurvature power spectrum. 

When H in f/27t > F a , the PQ symmetry restores during the inflation* 7 and in that 
case CDM axions are produced from the system of axionic string-wall system as well as 
the coherent oscillation with an initial alignment angle 9 = 0(1). This excludes H\ n i > 
10 n GeV for H ini /2n > F a in Fig. 2 (For details, we refer to e.g. Refs. [72, 73]). 

These constraints lead that even if 9 can be small by chance, large F a > 10 11 GeV is 
not allowed unless the inflation scale is low H in f < 10 11 GeV. 



6 Conclusion 

We studied constraints on non-Gaussianity in a mixture of adiabatic and uncorrelated 
CDM isocurvature perturbations. We adopted the optimal bispectrum estimator for the 
non-linearity parameters, and a fast method for non-Gaussian CMB simulation and the 
optimal weighting of observed maps are integrated in our analysis. Using the WMAP 

7 year data of CMB temperature maps with template marginalization of the Galactic 
foregrounds, we obtained a constraint <y 2 f^i°^ = 40 ± 62 at la level for the scale in- 
variant isocurvature power spectrum when the adiabatic perturbations are assumed to be 
Gaussian. Under the same setup, we also obtained a joint constraint on the non-linearity 
parameters (/nl, « 2 /nl°' > ) = (38 ± 24,-8 ± 72). The constraints weakly depend on the 
fiducial value of the isocurvature spectral index. Effects of the Galactic foregrounds at 
large angular scales are not severe. We found no statistically significant deviation from 
Gaussianity and the current WMAP observation of CMB is consistent with Gaussianity 
even when we include this type of isocurvature perturbations. We applied our results to 
the axion model. 

Since the CDM and baryon isocurvature modes affect the CMB anisotropy in the same 
way except for the overall amplitude, our constraint can be translated into the baryon 

isocurvature perturbations. This can be easily done by substituting (j^j a 2 f^°^ f° r 

# 7 Here we assumed the reheating temperature T rc h is below H in f . If T rc h > H in f , the PQ symmetry can 
restore for smaller Hi n f after inflation. 
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spectrum. See text for details. 
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a 2 /j~j L , where Qj, and f2 c are the density parameters of baryon and CDM, respectively. 
Although we restrict ourselves to uncorrelated CDM isocurvature perturbations in this 
paper, our method can be generalized to the cases of non-Gaussianity in isocurvature 
perturbations correlated with adiabatic ones and other types such as neutrino ones. These 
will be studied in a forthcoming paper [52]. 
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A Equivalence between {S^m) f B L =i an d (^prim^prirn)o 

In this appendix, we show that the normalization factor (Sp Tim ) fB =1 and the covariance 

(S , ^ rim 5' ] f rim )o is equivalent. 
From Eq. (15), we obtain 

(^prim)/^ L =l = q ^hhl3^mirn2m3 (^hmi^l2m2^l3m 3 ) f^ h =lJ (29) 

{lm} 

where we used the fact that (a/ m ) m _■, = 0. As discussed in Sec 3.3, we take di m = 

1 ' J NL x 

y^[C l ]lm,l'm' a l'm', SO that 
I'm' 



( a hmi a l 2 m 2 a hm3 ) /^ L =l 



(30) 

/j [C 1 ]l 2 m 2 ,l' 2 m' 2 [C %3m 3) l' 3 m' 3 (%m[ a l' 2 m' 2 a l' 3 m' 3 ) ^=1 

{I'm'} 

P ^imi^m'JC l \hm 2 ,V 2 m' 2 [C ,l' 3 m' 3 G m^m!,^^ 

{I'm'} 
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Combining Eqs. (29) and (31), we obtain 

(Sprim) fZ L =l = q ^ ^WsSmsP ^hmu^m'^C l ]l 2 m 2 ,l' 2 m' 2 [C Ibma^ra^^^^l^ 
{Iml'm'} 

(31) 

On the other hand, 



(^im'Sprini)o ~~ gg ^hhh^ l 2 l' 3 ^mim 2 m3^m[m' 2 m' 3 (^hmi^l 2 m 2 ^l 3 m3^l[m[^l 2 m' 2 ^l' 3 m' 3 )o (32) 

{iml'm'} 

— 6(az imi a^ m / i a^ m ^a z */ jm ^)o(ai 2m2 ai 3m3 )o + 9(ai imi a^ m / i )o(a; 2m2 a/ 3m3 )o(a^ m / ! a^ m / ! )o 
By taking the Wick's contraction, the terms in the square bracket are reduced into 

\C ]i 1 m 1 ,i' 1 m' 1 [C~ 1 ]i 2 m 2 ,i' 2 m' 2 [C~ ]i 3 m 3 j' 3 m 3 + (5 perms), (33) 
where we used the fact that (az m az' m ')o = Ci m y m i. Thus, we obtain 

( ^rim ^im)o = g ^ ^hhh ^fz 2 Z 3 ft«"3 ^m[m' 2 m 3 \P ^hrrii^m^ [C l \l 2 m 2 ,l' 2 m' 2 [C l ]l 3 m3,l' 3 m' 3 - 
{Iml'm'} 

(34) 

From Eqs. (31) and (34), we can finally see that (S£ rim ) fB _ =1 = (S^. ira S^ im ) . 



B Estimation of normalization 

Here we show how we can estimate the normalization factor (S^ m ) f B _-, from simulation. 

\ pnm/j NL — i 

First we divide the CMB anisotropy af m into Gaussian and non-Gaussian parts 

a hn = a G,lm + /NL a NG,Zm' (^5) 

where <2q /m and a^ G ; m respectively are the Gaussian and non-Gaussian contributions, 
which are given by 

4m = J drr 2 af(r) J dfY* m {r)X^{f), (36) 

o£g M = J drr 2 a?(r) j ^(f)(X^(f) 2 - <X£(f) 2 ». (37) 

Since a/ m is a linear combination of a/ m , we can also divide the weighted map into its 
Gaussian and non-Gaussian parts, af m = aQ lm + /nl^ng,;™- Then (S^ lim ) /s =i in Eq. (29) 
can be rewritten as 

(Sprim) f£ L =l = q bf ll2ls G^m 2m3 ((aG,hmi + ®NG,hmi ) (^G,Z 2 m 2 + ^NG,l 2 m 2 ) (^G,Z 3 m 3 + ^NG,Z 3 m 3 ) ) 

{lm} 

= k ^ 2 Z 3 ^mi 2 m 2 m 3 (°NG,/imi ^G,Z 2 m 2 ^G,Z 3 m 3 ) + (^NG,/ 2 m 2 ^NG,/ 2 m 2 ^NG,i 3 m 3 )) fB =1 > 

{Zm} 



(38) 
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where the last equality follows from the Wick's theorem. As long as non-Gaussianity is not 
too large, the second term in the last line can be omitted. Although we need many MC 
samples to evaluate the term (a§ G i imi aQ l2m2 aQ t ) , we checked that O(10) MC samples 
are enough for this purpose. In addition, since the terms in the last line of Eq. (38) consists 
of only anisotropy af m from a single perturbation mode X B , we can estimate (Sp rim ) /s =1 
without assuming a fiducial value for the fraction of the isocurvature contribution a. 
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